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ABSTRACT 

In the full-orbit particle simulations of energetic particle transport in plasmas, the plasma turbu- 
lence is typically described as a homogeneous superposition of linear Fourier modes. The turbulence 
evolution is, however, typically a nonlinear process, and, particularly in the heliospheric context, the 
solar wind plasma is inhomogeneous due to the transient structures, as observed by remote and in- 
situ measurements. In this work, we study the effects of the inhomogeneitics on energetic particle 
transport by using spatially distributed, superposed turbulence envelopes. We find that the cross-field 
transport is significantly reduced, when compared to the results obtained with homogeneous turbu- 
lence. The reduction can reach an order of magnitude when the enveloping breaks the wave phase 
coherence along the mean magnetic field direction. 
Subject headings: cosmic rays - diffusion - turbulence 
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1. INTRODUCTION 

When studying solar energetic particle (SEP) events, 
it is crucial to uncover the propagation of the particles 
in the interplanetary medium. Energetic particles related 
to solar eruptions are observed at large range of latitude s 
and longitudes fe.g.. iDalla et al.ll2003t iLiu et al.l 1201 If ) ■ 
and without understanding the propagation of particles 
across the mean magnetic field it is difficult to determine 
the extent of the acceleration region, and its relation to 
the physics of the eruption. 

The modelling of cosmic ray transport in turbulent 
plasmas i s typically ba sed on the diffusion-convection 
equation (Parker 1965), and the quasilinear approach 
incorporating pitch angle diffusion description, with 
the cross-fie ld propagatio n described as field line ran- 
dom walk (Jokipll 119661 ). The modelling of pitch 
angle e voluti o n has evolved to include adiab atic fo- 
cusing (lEarll I1976D. adiab atic deceleration (|Ruffolol 



Il995t iKocharov et al. |"P)9M and structured solar wind 
(|Kocharov et alj 12009) . The application of these mod- 
els to obtaining the SEP injection close to the Sun from 
the observations at 1 AU mostly concentrated on the 
particle propagation along the m ean magnetic field (e.g. 
lLaitinen et al.ll2000t [Drogc 2003) . Only recently has the 
cross-fiel d propagation been included into the S EP stud- 
ies (e .g. IZhang et al.l [200l iDroge et al.l [20Tot iHe et al.1 

[MB. 

A wide range of values for the perpendicular diffusion 
coefficient have been report ed. Studies of the t ransport 
of galactic cos mic rays (e.g. iBurger et al.ll2000T ) and Jo- 
vian electrons (jZhang et al.ll2007l and references therein) 
give estimates of order k±/k^ = 0.01. On the other hand, 
IZhang et~a l. (2003) report values of = 0.25 in re- 

lation to protons o riginating from a solar event, while 
iDwyer et al.l (|1997[ ) find the ratio to be of order of unity 
for energetic particles observed in association with coro- 
tating interaction regions. 

The values of k±/kh obtained from observations are 
related to different regions in the heliosphere, and possi- 
bly to different properties of the turbulence, which poses 
a problem when comparing with theoretical work. In re- 



cent years, the use of full-orbit simulations of energetic 
particles in studying particle propagation has gained at- 
tention. The benefit of this approach is that the param- 
eters defining the turbulence properties can be explic- 
itly prescribed. In addition, for the energetic particle 
propagation no a priori assumptions, such as the appli- 
cability of the diffusion des cription, are required. Th e 
approach has been used by Giacalo ne fe Jokipiil |l999), 
who obtained k±/k» — 0.01 — 0.03, with the perpen- 
dicular diffusion coefficient con tr adicting th e quasi linear 
theory results of I Jokipiil (|1966Q . iQ.in et al.1 (|2002D used 
the method to discover subdiffusion phase in the particle 
propagation perpendicular to the mean magnetic field, 
followed by subsequent recovery of diffusion. This phe- 
nomenon has subsequently been modelled as interplay 
between parallel and p erpendicular propagation effects 
(|Matthaeus et all 120031) . Thus, the full-orbit particle 
simulations have proved a valuable tool when studying 
particle transport. 

The full-orbit particle simulation results may depend, 
however, on the way the turbulent electric and mag- 
netic fields are described. A popular choice is to su- 
perp ose Fourier modes on a constant magnetic field 
fe.g.lGiacalone fc JokipiiirT^IQid^qollQin et alM M 
iZimbardo et al.l 120061: iRuffolo et al.l 12008) or m ore re- 
cently, on a Parker spiral field (|Tautz et al. 120111 ) . In this 
approach, however, the modes remain linear throughout 
the simulation domain, while the heliospheric plasma is 
known to evolve non-l inearly in addition to linear mode 
propagation (see, e.g-. lTu fc Marsch1ll995l for a review). 
Thus, the Fourier-mode description is not consistent with 
the behaviour of plasma turbulence, and using the ap- 
proach may result in inaccurate particle diffusion coeffi- 
cients. Some work exists on simulating particles in tur- 
bulent fields obtained by using MHD simulations, where 
the nonlinear interact ions are properly addressed (e.g., 
iBeresnvak et al.l 1201 1| ). However, such studies are lim- 
ited by computational requirements to a small range of 
scales and to cartesian geometries with constant back- 
ground fields. 

The assumption that the plasma turbulence is spatially 
homogeneous, while useful when comparing the simu- 
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lations to theoretical models, can also pose a problem 
when trying to understand SEP observations. The helio- 
spheric magnetic field expands as a Park er spiral, and 
may have a complicated stru cture (e.g., iParkerl Il958t 
iFiskl 119961 : iSternal et all 120 111 ). The solar wind is per- 
meated by transient structure s, as shown by in-situ a nd 
remote-sensing observations (jRouillard et al.l l2010al lbh. 
thus making it a spatio-temporally varying medium, in- 
stead of a steady solar wind flow. Thus, for SEP analy- 
sis there is a need to understand how such a structured 
plasma environment influences energetic particle propa- 
gation. 

In this work, we study the effect of structures in the 
turbulence on energetic particle propagation by means 
of full-orbit energetic particle simulations. We describe 
the turbulence as a superposition of localised, randomly 
distributed envelopes that contain a sum of linear waves. 
For simplicity, the envelopes are not limited in the direc- 
tion perpendicular to the mean field, but the variation 
takes place only along the mean field direction. We in- 
troduce two different models to study the effects of the 
structures. In the modulated wavefield model, only the 
amplitude of the turbulence is modulated, with the indi- 
vidual Fourier modes remaining coherent throughout the 
simulation domain. In the random envelope model, the 
phases, wave normal directions and polarisations of the 
modes are different in different envelopes. This is done to 
mimic the non-linear evolution of turbulence being con- 
verted throughout the heliosphere. Of these two models, 
the random envelope model captures better the influence 
of nonlinear interactions present in the heliospheric tur- 
bulence. Comparing the results of the random envelope 
model with those of the modulated wavefield model al- 
lows us to distinguish between the effects of amplitude 
modulation on one hand, and the phase, wave normal 
and polarisation changes on the other, as both effects 
are present in the random envelope model. 

The turbulence is superposed on a constant back- 
ground magnetic field. While such a background field 
limits the applicability of the model results to the he- 
liospheric environment, it is useful as a first approach 
in isolating the energetic particle transport effects due to 
the structured turbulence. The effects of a spatially non- 
uniform heliospheric magnetic field, and of fully three- 
dimensional structures, will be the subject of future stud- 
ies. 

In Section 2 we introduce our description of the struc- 
tured turbulence, and our approach to obtain the en- 
ergetic particle diffusion coefficients. In Section 3, we 
present the turbulence characteristics, and the full-orbit 
particle simulations and the derived diffusion coefficients, 
and study their variation with envelope characteristics, 
which can be considered as a measure for the structure 
size in the turbulent plasmas. We discuss the implica- 
tions of the results and draw our conclusions in Section 4. 

2. MODEL 

2.1. Turbulence Model 

The turbulence model consists of envelopes containing 
fluctuating magnetic field given by a sum of infinite plane 
waves. The magnetic field is defined as 

B = B a e z + 8B(x,y,z), 



where the background magnetic field is uniform and con- 
stant, with Bo = 5 nT, consistent with the magnetic field 
at 1 AU. The fluctuating fie ld <5B is based on the homo- 
geneous turbulence model of lGiacalone fc Jokipiil (|1999f ). 
and is defined as a sum of N modes, 



SB(x,y,z) 



N 

E 



A(k n )i n exp{i(k n z' n + /?„)}. (1) 



Here is the polarisation vector, 



cos a„x„ + i sin a n y n . 



The unprimed coordinate system has the background 
magnetic field along the z-axis, and the primed coordi- 
nate system is obtained from the unprimed through ro- 
tation characterised by the spherical coordinate system 
angles n and 4> n o f the w ave vector k„ , as defined in 
iGiacalone fc Jokipil (|1999j ). The fluctuation amplitude 
A(k n ) is given by a power law spectrum 



A 2 (k n ) =B\ 



G(k n ) 



AV n 



1 + {k n L c y 



where B\ is the variance of the magnetic field, L c the 
spectrum's turnover scale, 7 the spectral index, and AV n 
specifies the volume element in fc-space that the discrete 
mode k n represents. For the turnover scale, also called 
correlation length, we use L c = 2.15 r©, where r is the 
solar radius. 

In our simulations, we use the composite spectrum 
model, which is composed o f a sum of slab and two- 
dimensional (2D) turbulence (|Grav et al.lll996l ). For the 
slab component, k„ is directed along the background 
field, which is parallel to the z-axis, and the polarisa- 
tion vector lies in the xy-plane, with a random polari- 
sation angle a n . For the 2D-component, k„ lies in the 
xy-plane. The polarisation vector is also in the a;?/-plane, 
with a n = 7r/2, perpendicular to k„, to satisfy V-B = 0. 
The azimuthal angle <f) n and the random phase /3 n are 
chosen from a uniform random distribution. 

The slab and 2D spectra have the same turnover scale, 
and the wavemodes are generated for the same wavenum- 
bers k n . The wavenumbers are logarithmically spaced 
between 27r/lAU and 27r/10 _4 AU. The spectral indices 
for the slab and 2D components are 5/3 and 8/3, respec- 
tively, and the factor AV n is Ak n for the slab component 
and 2irk n Ak n for the 2D component. The total turbu- 
lence energy is divided between the slab and 2D compo - 
nent with ratio 20%:80%, following iBieber et al.1 (|1996l ). 

2.2. Turbulence Envelope Model 

We envelope the infinite plane waves, given by Eq. (fTJ 
into N p wave packets of length L p , randomly distributed 
in the simulation box. The envelopes are of cosine shape, 
given by 



2UC0 = 



1 



1 — cos 2n 



L, 



(2) 



e{z - zi} - e{z- (zi + Lp)} 



with Zi the starting point of the envelope, and the Heav- 
iside functions Q{z} limit each envelope to the range 
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Zi < z < Zi + hp. The magnetic field in these en- 
velopes remains divergence-free, as the envelope ampli- 
tude changes only in direction perpendicular to the polar- 
isation vectors £„. The overall turbulent magnetic field 
is the sum of the magnetic field in the envelopes, 

N p 

5B(x,y,z)=J^QL i (z)5B i , (3) 
i=i 

where <5B.; is given by Eq. ([1]). The distribution and sum- 
ming of the envelopes is shown in the bottom panels of 
Fig. [T] where the grey-filled curves show the amplitudes 
of each envelope, and the thick curve their sum. 

The envelopes are spread randomly along the z- 
direction from z = to z — L max , where the length of 
the simulation volume along the z-axis, L rnax , is chosen 
large enough that the particles will remain in the volume 
for the duration of the simulation. As the variation of the 
amplitudes takes place only in the z-direction, there is no 
need to limit the volume in the ccy-plane. The number of 
the envelopes, N p , is determined by defining an envelope 
density p p = N p L p / L max , which describes the filling of 
the simulation space with the turbulent envelopes. 

We consider two different models for <5B.; . In the wave- 
field modulation model, the set of angles /3„, 4> n and a n 
are the same for each envelope, and thus the enveloping 
only modulates the turbulent field described by Eq. ([1]). 
In the random envelope model, each envelope i has a 
unique set of angles /3„i, <p n i and a n i, mimicking the 
spatial evolution of the turbulence phases and wavenum- 
bers. 

2.3. Turbulence Energy Density 

To compare the turbulent fields and their effects on 
SEP propagation, we scale the fluctuation amplitudes so 
that the average energy density in the fluctuating field is 
independent of the enveloping parameters. As the energy 
density varies considerably in the enveloped cases, we 
calculate the energy density over a large spatial domain. 
The average energy density in a volume of cross-section 
S and height L max is given by 

U = -!-— / / SB 2 dz dS=^- (SB 2 ) , 

8^ S L max J J 8tt n 1 

where (SB 2 ) is the variance of the fluctuations over the 
volume. The energy density for the enveloped turbulence 
is then 

x r L max N p N p 

U env = — / V V KiKjSBi ■ SB, dz. (4) 

STTL max JO 

For the random envelope model, this integral can be es- 
timated analytically (see Appendix [A}. For the wavefield 
modulation model, the scaling factor is obtained through 
numerical integration. 

2.4. Energetic Particle Simulations 

In order to study charged particle propagation in tur- 
bulent magnetic fields, we integrate the fully relativistic 
equation of mot ion of energetic pr o tons using the sim- 
ulation code by iDalla fc Browning) (|2005l ). As the fluc- 
tuating field is magnetostatic, the electric field is zero. 



For the integratio n, we use the Bulirsh-Stoer method 
(|Press et alj|1993l) . The method uses adaptive timestep- 
ping to control the accuracy by limiting the error between 
consecutive steps to a given tolerance. We use the toler- 
ance of 10 -9 , which ensures the conservation of energy 
to within 10 -6 in our simulations. The particles' initial 
velocity distribution is isotropic and monoenergetic, and 
their initial positions are selected randomly in a region 
of size 40 times the turbulence correlation scale, in or- 
der to minimise the effects of local magnetic field. For 
each set of turbulence and particle parameters, we run 
the simulations for 10 field realisations, with 2048 parti- 
cles for each simulation. In the turbulence realisations, 
we typically use N = 128 Fourier modes, with validation 
runs using N = 1024 modes to verify that the results are 
not altered by the number of modes. 
The diffusion coefficient is obtained from the definition 

^CC = 2t ' C^x,y,z, (5) 

(e.g. IGiacalone k. Jokipiilll 999ft . where the square of the 
displacement, A£ 2 , is averaged over the simulated parti- 
cles. The diffusion coefficients are calculated separately 
for each field realisation. In Figs. [3] O K par refers to k zz , 
whereas n perp is obtained as the mean of k xx and n yy . 

The diffusion coefficients experience typically super- 
and subdiffusive phases, after which they settle to con- 
stant values, which we use in this study. In practice, 
we continue the simulations to ~ 100 parallel diffusion 
times in each run. In addition to this, we require that the 
particles have spread over a large number of envelopes. 
For this, we require that the FWHM of the distribution, 
as calculated from the analytic solution to the diffusion 
equation (e.g. IGiacalone fc Jokipiil | l~9 99), 4y/K zz tln2 is 
at least 40L p . This is done in order to reduce the statis- 
tical errors caused by the local differences between dif- 
ferent field realisations. 

We have reproduc e d som e of the simulations done by 
IGiacalone fc JokipiH (|1999ft . and found good agreement 
with their example run presented their Appendix A. We 
found, however, that the agreement between their ex- 
periments 1 and 3 is somewhat dependent on the se- 
lection of simulation domain size and fitting parameters 
when using the experiment 1. For this reason, the results 
presented in this study may differ somewhat from their 
results, as they use the experiment 1 for their results, 
whereas our method corresponds to their experiment 3. 

The diffusion of the particles can be caused by fieldline 
meandering, and by the drifts of the particles relative 
to the fieldlines. However, these processes are strongly 
coupled. Even a small drift of a particle from its original 
fieldline can result in large deviation from the original 
fieldline, depending on the local structure of the magnetic 
field. For this reason, separating the contributions of the 
two effects is not trivial. In this study, we do not attempt 
to make this distiction, but consider the diffusion as a 
compound effect of the two processes. 

3. RESULTS 

In this work, we study the effect of turbulent structures 
on energetic particle transport by using the turbulence 
model described in Section [2] varying the enveloping pa- 
rameters. For the envelope density p p we have chosen 
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Fig. 1. — Top panel: The sum of envelope amplitudes for pa- 
rameters p p = 16 and L p = 4L C (solid curve), and p v = 4 and 
L p = 1L C (dotted curve), in the modulated wavefield model. The 
horizontal dashed line represents the homogeneous model. Here, 
the energy density scaling factor 4/p 2 has been used. Bottom pan- 
els: The curves filled with greyscale colours show the individual 
envelopes, with L p = L c . The thick curve represents the sum of 
the envelopes. 

to use values 1, 4 and 16, to represent tenuous to dense 
enveloping. The length of the envelope, L p , is chosen 
to have values L c , 4L C and 16 L c . In this section, we 
first study how the enveloping affects the turbulence, and 
then how the particle propagation is affected. 

3.1. Enveloped Turbulent Field 

To understand how the enveloping affects the turbu- 
lence, it is useful to compare the behaviour of the en- 
velope amplitudes to the homogeneous turbulence. In 
the top panel of Fig. Q] we plot the sum of the enve- 
lope amplitudes for the wavefield modulation model, for 
two combinations of p p and L p . As can be seen, the 
high-density enveloping, with p p = 16, results in smaller 
deviations from the homogeneous model (the horizontal 
dashed curve) than the intermediate density (p p = 4). 
Thus, the envelope density describes the strength of the 
deviation of the mean amplitude from a constant value. 
The tenuous enveloping, with p p — 1 (bottom right panel 
of Fig. [1} , implying an average of one envelope at a given 
position, results in regions with no turbulence at all, and 
regions with overlapping envelopes. While such a sit- 
uation is not necessarily realistic, it is taken as an ex- 
treme of the parametrisation, and can in later studies 
be more realistically supplemented with some small un- 
derlying homogeneous turbulence. At the large density 
limit, the wavefield modulation approaches the homoge- 
neous turbulence description. 

We can also visualise how the field line wandering is 
affected by the enveloping. Fig. [21 shows the fieldlines 
for homogeneous turbulence (top panel), the modulated 
wave field (middle panel) and the random envelope model 
(bottom panel), for parameters L p = 4L C and p p = 4. 
The field lines are initiated from an area of 10L C x 10L C 
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Fig. 2. — The field lines in the homogeneous (top panel), mod- 
ulated wavefield (middle panel) and the random envelope turbu- 
lence model (bottom panel). In the latter two models, the envelope 
length is L p = AL C and density p p = 4. 
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Fig. 3. — Power spectrum of the turbulence using the homo- 
geneous (top panel), modulated wavefield (middle panel) and the 
random envelope turbulence model (bottom panel). In the latter 
two models, the envelope length is L p = AL C and density p p = 4. 

at z = 0. The most striking difference can be seen in the 
reduced spreading of the field lines in the case of ran- 
dom enveloping (bottom panel of Fig. [2]). There are some 
qualitative differences also between the homogeneous and 
modulated wave field models (the top and middle pan- 
els). However, they cannot be quantified on the basis 
of one realisation, and we will study their difference by 
using the energetic particle simulations in Section 13.21 

It is also important to determine how the turbulence 
spectrum is affected by enveloping. This can be seen in 
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Fig. [3] where we present a 2-dimensional Fourier trans- 
form of the turbulence, for the homogeneous (top panel), 
the modulated wave field (middle panel) and the random 
envelope turbulence model (bottom panel) . We show the 
spectrum for waves with the k along the mean magnetic 
field (solid line), along an axis perpendicular to it (dashed 
line), and the along direction 45° away from the mean 
magnetic field (dotted line). It can be seen that for the 
homogeneous turbulence the wave power resides in the 
parallel and perpendicular components only, whereas the 
enveloping effectively isotropises the turbulence at scales 
larger than the envelope length (L p = 4L C ). In the ran- 
dom envelope model, the perpendicular component of the 
spectrum (the dashed line in bottom panel of Fig. [3]) is 
reduced compared to the two other models. This may 
have an influence on energetic particle transport. 

It should be noted that in this Section we have pre- 
sented the behaviour of the amplitudes, fieldlines and 
spectra for only one realisation, and one region. Thus, 
the presented behaviour can be only be considered as 
qualitative, as localised structures are apparent even 
in ho mogeneous turbulence models (e.g. fChuvcha i et al.l 
l2007h . 

3.2. Energetic Particle Propagation 

To study the influence of enveloping the turbulence on 
energetic particle propagation, we have simulated sev- 
eral different scenarios, varying the particle energy and 
the turbulence properties. In Fig. [4] we show the par- 
allel diffusion coefficients obtained from these runs, for 
protons of energies 1 MeV, 10 MeV and 100 MeV. The 
points represent the mean, and the error-limits the stan- 
dard deviation of the diffusion coefficients, obtained from 
10 turbulence realisations. The top panel shows the mod- 
ulated wave field results and the bottom one those of the 
random envelope model. 

In both panels of Fig. |H the solid black line depicts 
the diffusion coefficient obtained from runs with homoge- 
neous turbulence, given by Eg . (Ill) . wh i ch are consistent 
with the results of lGiacalone &: Jokipiil ([1999). The non- 
solid lines show the parallel diffusion coefficient for en- 
veloped turbulence, with intermediate envelope density 
(pp = 4), and three envelope lengths. The enveloping 
reduces the diffusion coefficients somewhat for both the 
modulated wavefield model up to factor of 2 compared to 
the homogeneous turbulence model (top panel) and the 
random envelope model up to factor of 3 (bottom panel), 
and has a dependence on the envelope length. 

For the perpendicular diffusion coefficient, the effect of 
the enveloping is more varied. In Fig. [SJ we show the 
variation of the perpendicular diffusion coefficient as a 
function of energy for the same parameters as in Fig. [5j 
In the modulated wavefield model, the perpendicular dif- 
fusion coefficient is reduced by factor of two compared 
to the homogeneous turbulence result, and is indepen- 
dent of the envelope length. For the random envelope 
model, the reduction is strongly dependent on the enve- 
lope length, and the reduction in the diffusion coefficient 
reaches factor of 8 for short envelopes. 

The effect of the enveloping parameters can be seen 
in more detail in Fig. [51 where we show the perpendic- 
ular diffusion coefficient for three envelope densities, as 
a function of the envelope length. For the modulated 
wavefield model (top panel), there is no variation of the 
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Fig. 4. — Parallel diffusion coefficient as a function of particle 
energy, for the modulated wavefield model (top panel) and for the 
random envelope model (bottom panel), for p v = 4 and four enve- 
lope lengths. The solid line represents the result for homogeneous 
turbulence. The turbulence amplitude is B^/Bq = 1. 
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Fig. 5. — Perpendicular diffusion coefficient as a function of en- 
ergy, for the modulated wavefield model (top panel) and for the 
random envelope model (bottom panel). The turbulence ampli- 
tude is B\ /Bq = 1. 

perpendicular diffusion with the envelope length, within 
statistical errors of our study. There is a clear depen- 
dence on the envelope density, however, with a reduction 
of the diffusion coefficient up to factor 3 compared to 
the homogeneous turbulence result, when the turbulence 
enveloping is tenuous (p p = 1 in top panel of Fig. [BJ). 
Higher envelope densities result in the diffusion coeffi- 
cient approaching the homogeneous turbulence values 
When the phases, propagation directions and polari- 
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Fig. 6. — Perpendicular diffusion coefficient as a function of the 
envelope length, for the modulated wavefield model (top panel) and 
for the random envelope model (bottom panel). The turbulence 
amplitude is Bf/B^ = 1, and the particle energy 10 MeV. 

sations are different from one envelope to another, the 
cross-field particle propagation depends considerably on 
both enveloping parameters (bottom panel of Fig. |6]). 
For a high-density, long-envelope packeting, the limit of 
homogeneous turbulence is again obtained. When the 
envelopes are shortened and the envelope density is re- 
duced, however, the reduction of the perpendicular dif- 
fusion coefficient is nearly an order of magnitude. 

4. DISCUSSION AND CONCLUSIONS 

In this work, we have studied the effect of structured 
turbulence, as opposed to homogeneous turbulence, on 
energetic particle transport. We find that in the struc- 
tured turbulence the particle propagation is significantly 
inhibited, as compared to the propagation in homoge- 
neous turbulence of similar mean energy density of the 
turbulent magnetic field. 

The parallel diffusion can be seen to be affected 
strongly only when the inhomogeneities are of small 
scales (the dotted curve in Fig. @|. The reduction of 
the coefficient has also a significant energy-dependence, 
with higher-energy particles being affected more strongly. 
This is may be related to the increasing resonant scales 
of the particles (with Larmor radius of 0.4 r Q for 100 
MeV protons in a 5 nT magnetic field). 

The effect of the structuring of the turbulence on 
the cross-field propagation depends strongly on how the 
enveloping is performed. In the modulated wave-field 
model, where only the amplitude of the turbulence is 
modified, the cross-field propagation is sensitive only to 
the envelope density. For short envelopes, it appears that 
the effect of reduced field line wandering in the regions 
with small-amplitude turbulence is not fully compen- 
sated by the increased wandering in the higher-amplitude 
regions, but overall the cross-field propagation is inhib- 
ited. With the increase of the envelope density, the vari- 
ation of the turbulence amplitudes approaches the homo- 



geneous description (see Fig. [T]), and the limit of cross- 
field propagation in the homogeneous turbulence is ap- 
proached with increasing envelope density, as can be seen 
in the top panel of Fig. [5] 

For the random envelope model, however, the effect of 
the enveloping to the cross-field propagation is more var- 
ied, and stronger, and both the length and the density 
affect the perpendicular diffusion coefficient. Thus, the 
loss of phase coherence along the mean field direction, 
which is the basic difference between the wavefield mod- 
ulation and random envelope models, affects significantly 
to the cross-field propagation of energetic particles. 

This effect can be understood qualitatively. Strong 
cross-field wandering of a field-line, as seen in the top 
panels of Fig. [2 can only take place due to persistent 
cross-field disturbance with little or no dependence on 
the coordinate along the background magnetic field. This 
is possible due to the 2D component of the turbulence: a 
2D wave has a constant phase along the background mag- 
netic field and the cross-field direction perpendicular to 
the fc-vector. In the homogeneous and modulated wave- 
field model such waves are coherent throughout the simu- 
lation region (subject to disturbances from other waves), 
resulting in strong wandering of some field lines. In the 
random envelope model, however, this coherence is bro- 
ken, and such wandering can remain coherent only when 
one envelope dominates the others. 

Such breaking of the behaviour has natural justifica- 
tions, particularly in the heliospheric context. As dis- 
cussed in the Section [TJ the heliospheric magnetic field, 
and thus the turbulence, has its origins at the Sun, and 
the heliospheric plasmas are known to be structured due 
to transient phenomena at the Sun. Thus, it does not 
not seem feasible that the evolving corona could produce 
structures with infinite phase coherence into the helio- 
spheric magnetic fields. 

The non-linear evolution of turbulence also favours the 
evolution of the turbulence in parallel direction. As 
waves propagate outward from the Sun, they interact 
nonlinearly, resulti ng in cascade of ene rgy between dif- 
ferent scales (e.g., ITu fc Marschl [T995h . Although this 
energy cascade is very anisotropic, towards small perpen- 
dicular scales, the parallel scales are also affected. This 
can be understood t hrough the concept o f crit ical bal- 
ance, introduced by iGoldreich fc Sridharl (|1995fl . where 
the turbulence separated by the critical scale parallel to 
the mean field experience differing interactions due to the 
evolution of the counterstreaming waves. This results 
in scale-dependent anisotropy between the perpendicu- 
lar and parallel scales of the turbulence, which has been 
obser ved in the solar wind turbulence (Horbury ct al. 
l200l) . 

Overall, our results suggest that inhomogeneities 
within heliospheric plasmas and the nonlinear behaviour 
of the turbulence may significantly influence energetic 
particle transport, This has significant implications on 
the studies of SEP origins and their propagation in the 
heliospheric magnetic fields. 
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APPENDIX 

CALCULATION OF THE MEAN TURBULENCE ENERGY DENSITY 

For the random envelope model, the energy density given by Eq. (|4]) can be estimated analytically. As the magnetic 
fields 5Bi within different envelopes are not correlated, we can consider the integral over terms with i ^ j negligible. 
With this approximation, the integral can be written in the form 

1 1 N ? fL ma * 

U env = — V / SB 2 %(z) 2 dz dS. 

Using SB 2 (x, y, z) ps {5b 2 ) on the scales where 21^ changes appreciably, and ignoring the envelopes that are only partly 
in the integration region, the integral results in 

TT — - Lp 7" j 

U env — q j ^infi 

o ^max 

where Ui n f is the energy density of the turbulence given by Eq. (|T|). Thus, the scaling factor for the enveloped 
turbulence energy density is 

Uinf 8 L max 8 

^7V P U e nv 3 NpLp 3/9p 

The magnetic field amplitude is scaled by the square root of this scaling factor. We have verified the validity of this 
estimate through numerical integration of SB 2 . 

For the wavefield modulation model, the amplitudes can simply be added together, as the wave phases and polari- 
sations are the same in all envelopes. At high envelope density, the sum approaches constant, thus the scaling factor 
for the energy density is simply 4/p 2 . On the limit of low envelope density, the envelopes do not overlap, resulting in 
the same expression for the energy density as in the random envelope model, and the scaling factor 8/(3p p ). However, 
as the envelope densities in this study are in the intermediate range between these asymptotes, we use numerical 
integration for determining the scaling factor in the wavefield modulation model. 
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